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SIMPLIFIED INTERPLANETARY GUIDANCE PROCEDURES USING 
ONBOARD OPTICAL MEASUREMENTS 


By Harold A. Hamer and Katherine G. Johnson 
Langley Research Center 

SUMMARY 


Simplified guidance procedures are developed which are based on preflight deter- 
mination of the characteristics of trajectories perturbed about the nominal. The results 
are devoted principally to planetary approach guidance; however, some considerations 
for midcourse guidance are included. The methods are studied for an Earth -to -Mars 
trajectory but would be applicable to Grand Tour and other types of missions. 

Generally requiring only a single onboard optical angular measurement, the approach 
procedure predicts guidance corrections for the control of periapsis radius as well as 
orbital plane orientation. An error analysis was performed with an assumed la error 
of 10 seconds of arc in the optical measurement and a la velocity -cut -off error of 
0.1 m/sec. The analysis showed that when the usual type of midcourse guidance is applied, 
the periapsis radius at Mars can be controlled to a la accuracy of about 20 km if the 
approach guidance is executed 1/2 day before periapsis passage. Also, the orbital plane 
orientation can be determined to a la accuracy of 0.05°. If the guidance is performed 
at the Martian sphere of influence (2.2 days before periapsis), the periapsis -radius error 
doubles because of the increased effects of measurement error and maneuvering error. 

The effect of using another type of midcourse -guidance procedure was also investi- 
gated. This procedure differs in that the guidance correction and ensuing velocity errors 
are assumed to be approximately normal to the flight path. The analysis showed that this 
technique resulted in higher approach -guidance accuracy in controlling periapsis radius. 
However, the fuel requirement was increased by an order of magnitude and the accuracy 
of controlling the plane orientation was decreased. 


INTRODUCTION 


In space missions, the navigation and guidance are normally accomplished by auto- 
matic procedures which employ Earth-based radar measurements. The inclusion of pro- 
cedures which utilize onboard measurements is a desirable feature for interplanetary 
flight because planet ephemeris errors can lead to unacceptable errors in the trajectory 
position relative to the planet when only Earth-based measurements are used. This posi- 



tion error is most significant during the phase of the mission when the spacecraft is 
within the sphere of influence of the planet. Consequently, this paper concentrates on the 
approach-guidance problem. 

Over the years, a number of studies have been made to develop onboard guidance 
procedures for controlling the approach to a celestial body. (For example, see refs. 1 
to 7.) In general, these methods require several types of measurements and the mea- 
surements must be repeated a number of times inasmuch as the guidance correction is 
ordinarily based on statistical filtering techniques. The methods may also require more 
than one guidance maneuver. 

The approach procedure presented in this paper requires, in general, only a single 
onboard star-to-body angular measurement to determine the guidance correction. This 
measurement, made at a preselected time, is used to determine the position deviation 
from the nominal in a particular direction. When used in conjunction with some approxi- 
mations derived from two-body theory, a knowledge of this deviation is sufficient for con- 
trolling the periapsis magnitude by using a single guidance maneuver. The method is 
applied to the approach phase of a Martian trajectory. An analysis is included which 
shows that for reasonable accuracy, the deviation must be in a specified direction with 
respect to the nominal trajectory. Also, the direction of the guidance -correction vector 
is always parallel to the nominal orbital plane and perpendicular to the nominal velocity 
vector, which is essentially optimum. The accuracy characteristics of the method are 
examined by means of a Monte Carlo error analysis. The analysis includes the effects 
of measurement error, thrust-cut -off error, ephemeris error, and the approximation 
error caused by two-body assumptions and nonlinearity effects. 

Results are also included which show that the orbital plane orientation can be readily 
determined from an onboard measurement. Although plane changes are small at Mars, 
they can become important in a Grand Tour mission. 

Because of less restrictive midcourse -guidance accuracy requirements. Earth- 
based measurements will ordinarily suffice for this phase in most interplanetary missions. 
Brief results, however, are presented which could apply to Grand Tour missions, wherein 
midcourse guidance might necessarily be based on measurements relative to other planets. 

SYMBOLS 

A angle formed at spacecraft by lines of sight to centers of celestial bodies 

B angle formed at Earth center by line to spacecraft and line to Moon or Sun 

center 



C angle formed at Moon or Sun center by line to spacecraft and line to Earth 

center 

D position deviation in certain direction 

E position ephemeris error in direction of Sun 

R radius of Mars 

r range to Mars center 

r e range to Earth center 

r em distance between Earth and Moon centers 

r es distance between Earth and Sun centers 

r m range to Moon center 

r p periapsis radius at Mars 

r s range to Sun center 

Tp time to nominal periapsis time (time of approach -guidance measurement) 

V spacecraft areocentric velocity (Mars centered) 

Vp periapsis velocity at Mars 

x,y,z position coordinates in Cartesian axis system in which X-axis is toward the 

vernal equinox, XY-plane is parallel to earth equatorial plane, and Z-axis 
is in direction of north celestial pole 

j3 angle between AVp C and V (see appendix) 

y flight -path angle 

Ai difference in inclination angle of actual and nominal orbital planes 
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Ar difference between actual and nominal range to Mars center 

Ar e difference between actual and nominal range to Earth center 

Arp difference between actual and nominal periapsis radius 

AV approach -guidance velocity required to correct periapsis radius 

AV e difference between actual and nominal spacecraft geocentric velocity 

AVp C guidance velocity required for orbital plane change 

Ad angle between desired deviation direction and measurement star direction 

5 angle between deviation vector and orbital plane 

e eccentricity of orbit 

9 angle between deviation vector and body center 

0 m angle between measurement star and body center 

X angle between AV and V 

ji product of universal gravitational constant and mass of Mars 

a standard deviation or root -mean-square value of error 

cf) angle between Vj and V 2 (see appendix) 

Subscripts: 

A,B,C angles in triangle formed by vehicle, Earth, and Moon (fig. 1) 
a actual value 

D position deviation 

i inclination angle 


1 

t 


I 


4 



n 

nominal value 

R 

radius of Mars 

r,e 

range to Earth center 

r,P 

periapsis radius at Mars 

a 

half -angle subtended by Mars 

1,2 

value immediately before and after guidance correction, respectively (see 
appendix) 

Notation: 



| J absolute value 

A bar over a symbol indicates a vector. 

MIDCOURSE GUIDANCE 

The main purpose of this paper is to present results of applying a simplified onboard 
procedure for interplanetary approach guidance. Such a procedure corrects for errors 
incurred at midcourse; hence, a brief discussion of midcourse guidance is included. Only 
onboard midcourse guidance is considered here, inasmuch as this type of guidance was 
assumed for the approach guidance analysis. 

Previous studies of interplanetary midcourse guidance with onboard measurements 
(for example, refs. 8 and 9) have generally applied to manned flight wherein numerous 
observations over an extensive portion of the trajectory and, perhaps, several midcourse 
corrections are required. Results discussed herein correspond to the simplified 
midcourse -guidance procedure developed in reference 10. This method, based on pre- 
flight determination of the characteristics of trajectories perturbed about the nominal, 
requires only one position fix, from which is derived the standard fixed -point -of -arrival 
guidance correction. The simplicity of the measurements makes the method adaptable 
for flights such as the Grand Tour mission, where the midcourse -position fix with respect 
to other planets may be required. As is shown, the errors at the aim point may be rela- 
tively large. However, these errors can be corrected if the aim point is selected near 
the sphere of influence of the target planet. 


Range Determination 


In reference 10, the position vector of a body relative to the vehicle (position fix) is 
determined by three star-to-body angles and an onboard range measurement determined 
also from angular measurements. In reference 10 the error in the range measurement 
is shown to be the dominant factor affecting the accuracy of the position fix and hence the 
accuracy of the midcourse guidance. Because of the importance of this measurement, the 
pertinent aspects of onboard optical determination of range are presented in this section. 

A basic method of determining the range to Earth from onboard midcourse optical 
measurements taken in interplanetary space is illustrated in figures 1 and 2. The quan- 
tities r em and r es are known. (The method illustrated in fig. 1 can be applied to any 
planet if the position of one of its moons is known with reasonable accuracy.) 


Moon 



Spacecraft 


Earth 


Figure 2.- Vehicle -Earth-Sun system. 


For a Mars trajectory at a time 5 days from the Earth, typical values of the angles A, 

B, and C in figure 2 are 56.63°, 122.91°, and 0.46°, respectively. In figure 1, the max- 
imum value of the angle A at 5 days is approximately 16°. 

Two angles of the triangle are measured and the sine law is applied to solve for the 
range r e . For example, from figure 1 


r e = 


sin C 
sin A Pem 


| 
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or 


_ sin(A + B) 
r e - sin A rem 

The angle A can be measured directly, and the angle B (angle C) may be determined 
by (1) measuring the declination and right ascension of the Earth (Moon) or (2) measuring 
two star-to-Earth (star -to -Moon) angles. (See ref. 10.) In determining the angles B 
and C in this manner, their nominal values must be known. 

Range -Determination Accuracy 

In the vehicle -Earth-Moon system (fig. 1) the range-error equation, assuming mea- 
surement of angle B, is obtained from the total differential 

. 9 *V j * 9r p 9r p , 

dr e = TA + SB dB + dr <=m 

If the measurements are made at a selected nominal time, the third term may be neglected 
because of the certainty of r em ; hence, 

, ar p A 9r p m 
dre = 1 f dA+^dB 


For the region in which r e varies linearly with A and B, 


Ar ( 


9£e 

aA 




or 


i* 

Ar e = - . l A - AA + r em (cos B cot A - sin B)AB 
sin a 

If AA and AB are assumed to be random uncorrelated errors in the angular measure- 
ments, then 

(°r, e) 2 = ^a) + fem( cos B cot A - sin B)ct b ] 2 (1) 

For 120° > B > 60°, 

a re ^a A 
r > e sin A A 
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This approximation can be made because of the combined effect of the relatively small 
value of r em and the trigonometric values of the angles. 

Similarly, the range -error equation, if angle C is measured, is 


2 2 \2 
( CT r,e) = ( r e cot A a A) + ( r e cot c CT c) 

For 120° > B > 60°, 

CT r , e ~ r e cot A a A 


( 2 ) 


f 
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Figure 3 shows the variation of range error a r>e with the angle B when star-to- 
Earth and star -to -Moon sightings are used, as calculated from equations (1) and (2). Each 
curve involves two angular measurements, A and either B or C. Angle B is deter- 
mined by the relative position of the Moon. Obviously, measurements should not be taken 
when the Moon lies near the line of sight to the Earth. The data in figure 3 correspond to 
a spacecraft position 5 days from Earth (r e ~ 1 460 000 km). The la errors in the angu- 
lar measurements were assumed to be 10 seconds of arc; that is, the error in angle A 
would be 10 seconds of arc but the error in angle B or C would be greater inasmuch 
as two star-to-body measurements are required to determine these angles. From an 
example in reference 10, an error of 10 seconds of arc in each of these two measurements 
produced component errors in B or C of about 14 seconds of arc or a total la error 
of about 20 seconds of arc. 



Figure 3.- Range -determination error from angular measurements in vehicle- 
Earth-Moon system with la errors of 10 seconds of arc. 
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Equations were determined for the vehicle -Earth -Sim system (fig. 2); however, 
because of the geometry involved, the range errors were prohibitive. For example, at 
5 days from the Earth, the lex range error a r e is approximately 18 000 to 20 000 km, 
depending on whether the Sun or the Earth is used for the star -to -body angular sightings. 
In sighting to the Sun, the error is attributed to the small value of the angle C, whereas 
in sighting to the Earth, the error is caused by the large distance to the Sun. The effect 
of the uncertainty in the value of the astronomical unit is negligible. 

Guidance -Velocity Requirement 

The data in figure 4 illustrate the midcourse velocity requirement. The results in 
the figure are the velocity deviations 5 days from Earth due to all combinations of errors 



Figure 4.- Perturbed-trajectory characteristics at 5 days from Earth. 

of ±5 m/sec applied in one or more of the injection-velocity components. The velocity 
data are shown as a function of the range deviation inasmuch as this quantity is the main 
contributor to the midcourse correction requirement. The magnitude of the midcourse 
velocity correction would be slightly higher than the velocity deviation to account for the 
distance lost (or gained) at the time of first midcourse. The data of figure 4 are used 
subsequently to estimate the magnitude of the AV error associated with a first mid- 
course guidance correction predicted from onboard measurements. 

APPROACH GUIDANCE 


The simplified onboard procedure for interplanetary approach guidance is applied 
to the approach phase of an Earth-to-Mars mission. The approach region is considered 


9 


as that part of the trajectory lying within the Martian sphere of influence. The guidance 
procedure was developed in reference 11 for control of lunar approach trajectories. 
Characteristics of the method are different when applied to interplanetary approach 
because of geometry, spacecraft speed, and distance to the target body. The approach 
guidance results are shown primarily for the case where the final midcourse maneuver 
is performed relatively close to Mars. 

General Procedure 

A schematic sketch of the approach -guidance geometry is presented in figure 5. A 
nominal trajectory with a periapsis altitude of 1000 km ^rp = 4388 km) was chosen. The 



Figure 5.- Approach-guidance geometry. 

objective is to control the periapsis magnitude of the perturbed (actual) trajectory to the 
nominal value. Provision can also be made to control the orientation of the orbital plane 
to the nominal as discussed in the appendix. The plane changes are shown to be relatively 
small at Mars; however, such changes can become important in a Grand Tour mission. 

The measured value of the deviation D can be employed to determine the magnitude 
of the guidance correction required to control both periapsis radius and orbital plane 
orientation. 

The equation given in figure 5 shows that D is determined from nominal and mea- 
sured values of the angle 9 and the range. Range can be determined optically by mea- 
suring the angle subtended by Mars or by other angular measurements. The direction 
selected for D has a large effect on the guidance accuracy. The effect of range mea- 
surement error on D can be minimized by choosing a deviation direction perpendicular 
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f 


} 


to the nominal range vector; that is, 0 n * 90°. For this value of 9 n , the deviation can 
be calculated with the relationship 

D « r n cos 6> a 

Thus, the error associated with the range measurement is eliminated. As is shown later, 
this direction is usually best. Under certain conditions, however, another direction could 
prove superior because of reduced scatter error or increased measurement sensitivity. 

Where the final midcourse maneuver is performed far from the Martian sphere of 
influence, a range measurement is required. 

In figure 5, the deviation D is shown in the direction of the measurement star; as 
stated, the direction of D is most important. In practice, if no suitable measurement 
star lies in the desired direction of D, another star may be substituted in the following 
manner. Choose a star which lies in the plane containing r n and the desired direction. 
By adding (or subtracting) A 9, the difference between the nominal values (see sketch), 
the value of 0 a can be obtained from 


where 9 m is the angle between the measurement star and body center. When the chosen 
star is not in this plane, the direction of the desired star must be rotated about r n into 
the plane of the measurement star. The conversion of the angular measurement and the 
calculation of D would then be based on this new direction. 


Measurement Desired direction 
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Guidance -Velocity Requirement 


is 


The guidance velocity required to change the periapsis radius to the nominal value 


AV = 


v|r 2 COS y cos(y + X) - r 2 n cos xj 


r p, n - 1,2 cos2 (r + X) 


r p,nj( r2 - r p,n) v2 sin 2 X + r 2 cos 2 (y + X) - r 2 n ^ 

r p, n " r2 cos 2 (y + X) 


1/2 


( 3 ) 


This equation was derived in reference 11 from two -body relations. The values derived 
with the alternate signs of the second term correspond to correcting to either side of 
Mars. The lesser magnitude of AV would ordinarily be chosen to assure that the 
proper direction around Mars was achieved. The AV required per meter of periapsis 
radius is shown in figure 6 as a function of distance and time to Mars. (To a good 



16 x 10 5 



Figure 6.- Approach- guidance velocity requirement in relation to range 

and time to periapsis. 
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approximation, AV is linearly related to periapsis radius.) The curve, computed by a 
simplified equation (ref. 11) by using nominal values and X = 90°, illustrates the rapid 
increase in the velocity requirement as Mars is approached. 

In equation (3) the flight -path angle y has the major influence on the magnitude 
of AV for a given guidance pointing angle X and guidance maneuver time. It will be 
shown in the following section that the value of y is highly correlated with the deviation 
taken in certain directions. A measurement to a star is used to determine D, which in 
turn predicts the AV requirement. 

Results With Close Midcourse Maneuver 

General considerations . - For the data presented in figures 7 to 13, an assumed 
final midcourse correction was performed close to Mars. The time selected was 10 days 
before arrival at Mars. A Monte Carlo procedure was employed to simulate midcourse 
guidance errors due to onboard measurement errors resulting from the type of onboard 
measurements previously discussed. The onboard measurement -error analysis of ref- 
erence 10 was used to determine the error distributions. (Although these distributions 
pertain to lunar missions, it was assumed that they would apply to interplanetary flight 
because of the similar measurements used.) The resulting guidance pointing error dis- 
tributions had la in-plane and la out-of-plane errors of 1.8° and 0.7°, respectively. 
The la error in AV was 2 m/sec. The effect of guidance maneuvering error was 
considered negligible. Except where noted, the midcourse errors were essentially along 
the flight path, which is the general case for most midcourse procedures. (For example, 
see ref. 11.) 

Examples of the variation of y with D are shown in figures 7 and 8. (The lines 
in these, as in all figures, were faired through the data points.) Each data point repre- 
sents the condition at Tp = 1 day on a perturbed trajectory resulting from the errors 
in the final midcourse correction. The deviation D was determined by using the 
relationship 

D = l (x a - x n ) + m(y a - y n ) + n(z a - z n ) (4) 

where l, m, and n are the known direction cosines of the deviation vector. Equa- 
tion (4) is equivalent to the equation in figure 5. In figure 7, the deviation lies in the 
nominal orbital plane and is perpendicular to the nominal radius vector, whereas in fig- 
ure 8, its direction has been changed 2° in the nominal orbital plane. 

The scatter of the data such as shown in figures 7 and 8 produces error in the guid- 
ance procedure. It is important that the amount of scatter be minimized by selecting the 
most effective direction for D. Comparison of figures 7 and 8 shows that this direction 

13 
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Figure 7-- Example of good correlation 
between 7 and D. £>n = 0°J 
T p = 1 day. 



Deviation, D, km 

Figure 8.- Example of poor correlation 
between 7 and D. B n = 0°; 

Tp = 1 day. 


should be approximately perpendicular to rp. The acceptable range for 0 n is approxi 
mately 89.5° to 90.5°. Varying the out -of -plane direction by 10° does not greatly affect 
scatter error; however, the scatter error tends to increase for angles much greater 
than 10°. 


The Sun can be used as the measurement star if it is reasonably close to the orbital 
plane of the spacecraft. The assumption that the Sun is an inertially fixed point in space 
contributes a maximum error in D of 1 km, which is negligible. Sighting on the Sun 
rather than on a star could prove far superior because of ease of acquisition. 

Figures 9 and 10 are presented mainly to indicate the amount of per iapsis -radius 
error caused by approximation (scatter) error. The distance between each point and the 
line represents the amount of error in determining periapsis radius for the corresponding 
trajectories. The la values in the figures show that somewhat higher accuracy is 
achieved with a deviation which is 10° out of the nominal orbital plane of the spacecraft. 
The values of the deviation D remain essentially unchanged for a given direction at any 
time within the sphere of influence of the planet. 

Uncorrected values of Arp in figures 9 and 10 are less than 200 km. These rela- 
tively small values are due to the fact that the errors in the final midcourse velocity cor- 
rection were generally in the direction of the trajectory motion relative to Mars. 




(< 


< 
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Deviation, D, km Deviation, D, km 

Figure 9 .- Accuracy of determining peri- Figure 10.- Accuracy of determining peri- 
apsis radius for 9 n = 90 °, 5 n = 0 °, apsis radius for 0 n = 90 °, S n = 10 °, 

and Tp = 1 day. and Ip = 1 day. 

To determine the guidance correction, only the variation of approach -guidance cor- 
rection AV with deviation D (fig. 11) is required. This variation is calculated with 
equations (3) and (4). In all cases, the AV vector is parallel to the nominal orbital 



Deviation, D, km 


Figure 11.- Guidance velocity as a function of deviation. 
9 n = 90°; S n = 10°; T p = 1 day. 
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plane and perpendicular to the nominal velocity vector, which is essentially optimum for 
the fuel requirements. Both the deviation and the AV values in figure 11 pertain to the 
same time, T p = 1 day. The guidance maneuver can be performed any time after the 
measurement through the use of the data given in figure 6. Note that a delay of 1/2 day 
doubles the fuel requirement. 

Effect of different midcourse procedure. - The preceding discussion has dealt with 
midcourse guidance corrections which produce an error pattern along the flight path. It 
is conceivable that midcourse errors in a perpendicular direction could be generated by 
some other midcourse guidance method. The effect of such an error pattern was studied. 
The data, in figures 12 and 13 represent approach conditions for perturbed trajectories in 
which the final midcourse errors were directed essentially perpendicular to V. The 
midcourse -guidance error distributions were assumed to be the same as those used to 
derive figures 7 to 11. 

Note in figure 12 that the range of Ar p has increased by an order of magnitude 
from that of figures 9 and 10 where the midcourse guidance errors were along V, but 
the scatter has been drastically reduced. In addition, the data in figure 12 show that a 
wide range of deviation directions is available for calculating D. The set of data shown 
for 5 n = 54.73° represents a measurement in the direction of the Sun. No increase in 

4xl0 3 



Figure 12.- Accuracy of determining periapsis radius with final midcourse 
errors approximately perpendicular to flight path. (Note staggered 
vertical scale.) 


16 
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Deviation, D, km 

Figure 13.- Effect of measurement time on deviation. 

9n = 90°; 6 n = 0°. 

the scatter error occurs because of this large out -of -plane angle. However, the sensi- 
tivity of D has been reduced, which, in turn, will increase the effect of measurement 
error on the guidance accuracy. 

In figure 13 it is seen that the sensitivity of D is not greatly affected by the mea- 
surement time. As will be shown, the ratio 9rp^8D is important in the error analysis. 
This ratio varies from about 1.06 to 0.90 when midcourse errors are perpendicular to V 
and the ratio remained constant at approximately 0.93 for all values of Tp when mid- 
course errors were along V (for example, figs. 9 and 10). 

Results Without Close Midcourse Maneuver 

For the data presented in figure 14, a first midcourse maneuver at a time 5 days 
from Earth was assumed to be the only guidance correction prior to approach guidance. 
As shown, this procedure leads to very large per iapsis -radius errors which must be cor- 
rected by an approach -guidance maneuver. Normally, the use of only one midcourse 
maneuver causes excessive fuel requirements. At the sphere of influence the lcr value 
of AV required for correcting the rp dispersions in figure 14 would be about 
120 m/sec. The fuel requirements could be lowered by applying the correction at an 
earlier time. (See fig. 6.) Even though the trajectory would lie outside the sphere of 
influence at this time, the two-body approximation used for equation (3) would still be 
adequate for the time period shown in figure 6. 
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Figure 14. - Deviations at Martian sphere of influence ^Tp = 2.2 days) 
resulting from first midcourse perturbations. 

A single midcourse correction might be justified for a Grand Tour mission where 
midcourse maneuvers between outer planets are based on measurements with respect 
to a planet other than the Earth. Here the highly accurate Earth-based radar measure- 
ments normally required for the more refined second midcourse correction would not be 
available. 

As in previous figures, each data point in figure 14 is the result of the trajectory 
having been perturbed at midcourse with a Monte Carlo procedure. The same midcourse 
error distributions as previously discussed were used, except that a lcr error of 1 m/sec 
in the AV magnitude was assumed. This value was based on data of figures 3 and 4. 

The minimum range -determination error from onboard measurements is shown to be 
about 250 km. This error is the main contributor to the midcourse guidance error and 
figure 4 shows that this amount of error would lead to a AV error of roughly 1 m/sec. 

The data in figure 14, with little scatter for the cases shown, indicate a wide range 
of usable deviation directions. The angle between the two planes referred to in the figure 
is about 55°. The small effect of direction on scatter error is attributed to the highly 
elongated shape of the position-error ellipsoid at the Martian sphere of influence. This 
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ellipsoid was determined by analyzing the errors in the perturbed trajectories at this 
point, which had resulted from the Monte Carlo midcourse -guidance errors. The 
la values of the three axes of the error ellipsoid are 265 000 km, 3880 km, and 550 km, 
which illustrates that most of the error is concentrated along the major axis. This axis 
lies about 9° out of the vehicle -Sun- Mars plane about midway between the 90° angle 
formed by the radius vectors to the Sim and Mars. The inherent advantage of choosing 
9 n ~ 90° is discussed in the following section. 

As aforementioned, the Sun can be assumed a fixed point in inertial space and used 
as the measurement star. This assumption causes a maximum error in the deviation D 
of only 10 km and, for 0 n = 90° (fig. 14), an error of about 10 km in determining rp. 

APPROACH-GUIDANCE ACCURACY 

In this section, the errors associated with the approach -guidance procedure are 
defined and analyzed and their effect on the accuracy of controlling periapsis radius is 
determined. The errors are summarized in table I for two times; the time Tp = 2.2 days 
is at the Martian sphere of influence. In addition, the effect of ephemeris error is shown 
in figure 15. The results in table I apply to a deviation direction taken perpendicular to 
the nominal range vector (d n = 90°). The guidance maneuver was assumed to be per- 
formed at Tp. 

TABLE I. - APPROACH-GUIDANCE ACCURACY CHARACTERISTICS 



a rjP , km, for - 

Type of error 

T d = 2.2 days 
(a) 

T D = 0. 5 day 

(a) i 

Measurement 

26.0 

6.2 



(29.5) 

( 6.0) 


Scatter 

17.2 

17.2 



( 2.0) 

( 2.0) 


Velocity cut-off 

17.0 

4.0 

Total 

35.5 

18.7 



(34.1) 

( 7.5) 



a Results in parentheses correspond to final midcourse -guidance errors 
in a direction approximately perpendicular to the flight path. 
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Error in astronomical unit, km 


Figure 15.- Effect of ephemeris error on approach-guidance accuracy. Tp = 1 day. 

The results in the table apply to the close midcourse case; that is, a final midcourse 
correction was employed several days prior to the approach guidance. The values are 
given for midcourse guidance errors both along and normal to V. The value for the 
scatter error was obtained from figures 10 and 12. This error is caused by the approxi- 
mation in predicting the guidance correction from a value of the deviation D. As shown, 
the scatter error is very small for the case where the midcourse guidance errors are 
normal to V. This, in turn, reduces considerably the total error for T p = 0. 5 day. The 
total error was determined by statistically combining the effects of the various types of 
errors. 

The la per iapsis -radius error due to maneuvering error (velocity cut-off) is 
shown for a typical la value of 0.1 m/sec. The periapsis error was determined by 
multiplying the reciprocal of the velocity requirement (fig. 6) by 0.1 X 10 The results 
in table I correspond to a guidance pointing angle X = 90°. By sacrificing the fuel 
requirement, which is relatively small (fig. 11), and by applying the approach AV vec- 
tor at X - 0°, the effect of the maneuvering error can be reduced significantly, provided 
that the maneuver pointing error is small. 

The effect of measurement error on the guidance accuracy can be determined from 
the following equations: 



20 



and 


<r r,p = W a D (6) 

where values for drp/dD are obtained from plots such as those given in figures 10 
and 12. Error equations for the onboard optical measurements were developed in ref- 
erence 11; equation (5) corresponds to uncorrelated errors in the measurement of 
range r and the angle 6. 

In equation (5), nominal values for r and 6 are used in calculating the error ctjj. 
For values of 0 n ~ 90°, it is seen that 


°D “ ra e 


(7) 


Equation (7), along with equation (6), was used to calculate the results shown in table I. 
When the midcourse errors were along V, the ratio of 8rp to 9D varied slightly with 
time Tp. (See fig. 13.) A typical la value of 10 seconds of arc was assumed for the 
angular measurement accuracy. 

For Q n * 90°, the effect of range -measurement error a r becomes important, as 
shown by equation (5). If range is measured by the half -angle subtended by Mars a, then 
according to reference 11 


Cr = 




+ r 


! 4 ) 


,1/2 


where aj^ is the uncertainty in the value for the radius of Mars and cr a is constant 
with time. As an example, at Tp = 0.5 day and with 0 n = 85°, the error in periapsis 
radius a r p due to measurement error is 28 km. (Typical values of 5 km and 10 sec- 
onds of arc, respectively, were assumed for cr^ and a a .) This error is due almost 
entirely to error in the range measurement and would be much higher at greater distances 
from Mars. Thus, if is much greater (or less) than 90°, some optical method other 
than the subtended -angle method would be highly desirable for measuring range. One 
point of import is that the la value of the incremental range Ar for the Monte Carlo 
perturbed trajectories may be well below the range -measurement accuracy. For exam- 
ple, the data for the close final midcourse maneuvers with guidance errors perpendicular 
to the flight path showed a la value of Ar less than 100 km at the Martian sphere of 
influence. In this case, then, regardless of the value used for 9 n , approximating range 
with its nominal value would be adequate. 
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Effect of Ephemeris Error 


The effect of Martian ephemeris error on the approach -guidance accuracy was 
examined (fig. 15). Ephemeris error affects the accuracy because the guidance measure- 
ments are referenced to a nominal trajectory which, in turn, is based on a certain location 
of Mars. This type of error was investigated only for position error in the radial direc- 
tion; error in the direction of the motion of Mars would have a negligible effect on r p . 

In figure 15, it is seen that the approach -guidance procedure compensates, to a large 
extent, for presence of ephemeris error. The data are shown for T p = 1 day but are 
representative of results for most times within the sphere of influence. The results, 
shown as the ratio of | Ar p j to the ephemeris error E, were obtained by using sev- 
eral trajectories with different nominal values of r p . Despite apparent differences in 
the results between the cases, the trend of the data indicates that increasing the error 
toward the Sun (which, in effect, moves Mars away from the nominal trajectory) does 
not increase the guidance error appreciably. The region shown in the figure for 
Ar p j < 10 km corresponds to the uncertainty (3a) in the Mars ephemeris which is in 
uhe neighborhood of 200 km (ref. 12). Adding the 10-km error, statistically, to the data 
of table I does not appreciably change the results in most cases. The excessive ephem- 
eris errors included in figure 15 are indicative of results which might be expected for 
other planets. 


Orbital -Plane Orientation Errors 

Results in the appendix show that the orbital -plane orientation errors can be con- 
trolled to high accuracy if deviations are measured in a certain direction. The acceptable 
ranges for these measurements, as well as for those used to control periapsis radius, are 
summarized in table n. As previously noted, the direction of the measurement star need 
not fall within these ranges. 

TABLE n.- APPROXIMATE ALLOWABLE RANGES FOR % AND S n 


Direction of 
final midcourse 
velocity vector 

Periapsis radius 
to be controlled 

Orbital plane orientation 
to be controlled 

9 n , deg 

6 n , deg 

9 n , deg 

6 n , deg 

Approximately 

89.5 

±10 

89.8 

88 to 90 

along flight 

to 


to 

and 

path 

90.5 


90.2 

-88 to -90 

Approximately 

70 

±60 

88 

89.8 to 90 

perpendicular 

to 


to 

and 

to flight path 

110 


92 

-89.8 to -90 
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CONCLUDING REMARKS 


A procedure has been developed for interplanetary approach guidance by using 
onboard optical measurements. The method is based on characteristics of perturbed 
trajectories. The procedure requires, in general, only a single star -to -body measure- 
ment to determine the position deviation from the nominal in a particular direction. An 
error analysis showed that the procedure was feasible for controlling the periapsis radius 
of Martian trajectories; periapsis radius can be controlled to a la accuracy of about 
35 km if the maneuver is performed at the Martian sphere of influence. The accuracy is 
substantially increased if the guidance is performed closer to Mars. 

Some preliminary results on midcourse guidance for interplanetary flight showed 
that midcourse procedures based on onboard measurements may be applicable to Grand 
Tour missions. 

The important results concerning the approach-guidance method are as follows: 

The method can be applied anywhere within the Martian sphere of influence which 
extends some 570 000 km. 

It has been shown that the guidance method compensates, to a large degree, for the 
effect of ephemeris error. 

The Sun may be assumed to be a fixed point in inertial space and used as the mea- 
surement star, with little degradation of guidance accuracy. 

Directions that may be selected for the deviation depend on the type of procedure 
used for final midcoin- se guidance. A direction lying perpendicular to the nominal radius 
vector and reasonably close to the orbital plane of the spacecraft, within ±10°, generally 
provides maximum accuracy for controlling periapsis radius. 

The type of final midcourse procedure has a pronounced effect on the approach 
accuracy, especially if the approach guidance is performed relatively close to Mars. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 13, 1972. 
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APPENDIX 


CONTROL OF ORBITAL PLANE ORIENTATION 

The approach -plane orientation errors at Mars are minor; however, if such errors 
at the outer planets are not corrected for Grand Tour missions, they could lead to extreme 
errors at the next planet. The purpose of this appendix is to show that these orbital -plane 
errors can be determined from onboard measurements. Results in figures 16 and 17 indi- 
cate that the error in the inclination angle can be predicted by a measurement of a devia- 
tion in a direction perpendicular to the nominal orbital plane. In figure 17, where the final 



Figure 16.- Variation of Ai with D for 9 n = 90° and S n = 90°. 

midcourse errors were perpendicular to the flight path, there is a noticeable scatter error 
with a lcr value of about 0.03°. As in the case shown in the main text for predicting r p 
from deviations taken near the orbital plane, the deviation direction for determining incli- 
nation angle may be rotated in a plane perpendicular to the orbital plane without affecting 
the scatter error as long as 9 n ~ 90°. The allowable rotation, however, is only about ±2°. 

The change in the location of the ascending node can be determined with the same 
degree of accuracy shown in figures 16 and 17, since in general 

A N = Aoo - sin _1 (tan cot i) 
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APPENDIX - Continued 



Figure 17.- Variation of Ai with D for final midcourse 
errors approximately perpendicular to flight path. 

0 n = 90 ° ? &n ~ 90 ° } Tp = 1 day • 


or 


dA N = (* ~ 


tan 2 ? c 


cot 2 i) 


- 1/2 


esc' 2 di 


where 

latitude of V at Martian sphere of influence 
Aoo longitude of V at Martian sphere of influence 

Ajj longitude of ascending node 

The change in Ajg- was found to be approximately 0.67 Ai in the present cases. 

In figure 16 the variation of Ai with D is not greatly affected by Tp. The scat- 
ter (approximation) error remains negligible at all values of Tp, except when Tp ~ 0. 
The guidance maneuver to correct the plane orientation (derived from the changes in the 
inclination and the ascending node) should therefore be made close to the planet to take 
advantage of the smaller effect of measurement error. As shown in figure 18, the 
guidance -velocity requirement is essentially independent of distance from the planet. 
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APPENDIX - Continued 



Time to nominal periapsis time, days 

Figure l8. - Guidance -velocity requirement for changing 
the orbital plane. 

The equation for AV pc was derived from the sketch 

Vl 



V 2 



These equations correspond to the case where the magnitude Vj = V 2 . (See ref. 1.) 

When the final midcourse errors are perpendicular to the flight path (fig. 17), 

Tp = 0.5 day is about the closest practical time at which to perform the plane-change 
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APPENDIX — Concluded 

guidance. At this time, the effect of the measurement error is roughly the same as the 
0.03° (la) effect of the scatter error. For example, 

°i = ~$b °b = °- 0073 "d 

where a^ at Tp = 0.5 day is approximately 6 km; hence, a^ ~ 0.043°. 
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